Theoretical prediction of multiferroicity in double perovskite Y 2 NiMnOe 
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We put forward double perovskites of the R^NiMnOg family (with R a rare-earth atom) as a 
new class of multiferroics on the basis of ab initio density functional calculations. We show that 
changing R from La to Y drives the ground-state from ferromagnetic to antiferromagnetic with 
Till spin patterns. This E*-type ordering breaks inversion symmetry and generates a ferroelectric 
polarization of few (iC/crn 2 . By analyzing a model Hamiltonian we understand the microscopic 
origin of this transition and show that an external electric field can be used to tune the transition, 
thus allowing electrical control of the magnetization. 

PACS numbers: 71.45.Gm, 71.10.Ca, 71.10.-w, 73.21.-b 
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Materials with simultaneous magnetic and ferroelec- 
tric ordering - multiferroics - are attracting enormous 
scientific interest due to their potential for applications 
in memory and data storage devices 0, [3] ■ Multiferroics 
with magnetically induced ferroelectric order are particu- 
larly interesting due to their strong magneto-electric cou- 
pling, which is required for an electric (magnetic) control 
of magnetic (electric) order parameter: a very desirable 
property from a technological point of view [3J. There- 
fore, the search for new multiferroics with the ferroelec- 
tric order driven by the magnetic order is currently a very 
active and important field of research. 

In this letter, we show, by using first-principles den- 
sity functional theory (DFT) calculations and a model 
Hamiltonian analysis, that the magnetic order in the dou- 
ble perovskite compounds R^NiMnOg (RNMO) changes 
from ferromagnetic for R=La and Sm to the E*-type for 
R=Y. The E*-type magnetic structure consists of T-T-|- 
l spin chains along the cubic perovskite-like directions 
(see Fig. \S§ , or equivalently, zig-zag ferromagnetic (FM) 
spin-chains antiferromagentically coupled in-plane to the 
neighboring zig-zag chain, with an out-of-plane FM cou- 
pling |4|]. E*-type magnetism breaks inversion symmetry 
and thus allows a ferroelectric polarization to occur. We 
will show that indeed Y2NiMnO@ is multiferroic with a 
electric polarization of few /j,C/cto 2 . We will show that 
vice versa the magnetic transition from ferro to E*-type 
can be tuned by an external electric field, thus allowing 
electric control of the ferromagnetic order parameter. 

The double perovskite La 2 NiMn0 6 (LNMO) is a fer- 
romagnetic insulator with a Curie temperature close to 
room temperature (T c ~ 280 K). The structure of LNMO 
changes from rhombohedral (R3) at high temperature, to 
monoclinic with P2i/n symmetry group at low temper- 
ature [1, 0]. DFT has been successful in understand- 
ing the ferromagnetic (FM) insulating behaviour and the 
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FIG. 1: (Color online) (a) Density of states in P2i/n sym- 
metry for FM-type magnetic structure in LNMO; (b) Density 
of states in P2i symmetry for E*-type magnetic structure in 
YNMO. Zero energy is set to be Fermi level. 



dielectric anomalies observed experimentally in LNMO 
0, Q ■ We work in the framework of density- functional 
theory, using the Vienna ab initio simulation package 0] , 
in which the Kohn-Sham equations are computed using 
the projector augmented wave method (PAW) and solved 
describing electronic exchange and correlation in the gen- 
eralized gradient approximation (PBE) [13, [HI ■ 

The spin resolved density of states (DOS) of LNMO 
is shown in Fig. [1] a) and they agree well with the one 
in Ref. 0]. The octahedral surrounding of Mn and Ni 
atoms split the Mn-d and Ni-d manifolds into the t 2g 
and e g levels. In the spin up channel the Ni-t2 9 and Ni- 
e g are found within the range of 2 eV below the Fermi 
level (Ep) mixed with O-p and the Mn-t2 9 states, the 
latter being localized between the Ni-t2 3 and the Ni-e s . 
The Mn-e 9 remain empty, separated by ^2.5 eV from the 
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Mn-t2 ff . In the down spin channel the Ni-t2 9 states are 
located between the O-p and the Fermi level, while the 
Ni-e s lie above the Fermi level along with the Mn-t2 9 and 
the Mn-e ff . This distribution of states corresponds to the 
nominal valences Ni 2+ (d 8 ) and Mn 4+ (d 3 ) : the mag- 
netic moments in our simulations are indeed found to be 
1.4 fib and 2.9 respectevely, close to the experimental 
values 
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Calculating the DFT total energies of the different 
magnetic structures on the relaxed crystal structure of 
LNMO and using a least square mean method to fit their 
values by an effective Heisenberg Hamiltonian (as in Ref. 
[Hi]), we extract the magnetic couplings Jy [l3 |. 
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1 f(Ni-Mn) denote the in-plane Ni-Ni second- nearest neigh- 
bor, in-plane Mn-Mn second-nearest neighbor, out-of-plane 
Ni-Mn nearest-neighbor and in-plane Ni-Mn nearest-neighbor 
exchange constants, respectively. The ratio (Q) between 
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is reported in the last column. 



Experimentally, the unit cell parameters of double per- 
ovskites R^NiMnOe with R = Sm, Y are unknown. Since 
the theoretical lattice parameters of LNMO obtained 
within GGA are very close to the experimental estimates 
(see also Ref. Q), we infer a good accuracy for the pre- 
diction of structural properties. Therefore, we determine 
the lattice parameters and the internal coordinates of the 
ions for R=Sm, Y by earring out the first-principles vol- 
ume and shape optimization for the monoclinic phases 
upon imposing a FM spin configuration. Our results 
may provide basis to refine the crystallographic struc- 
tures experimentally [l5j |. A general trend is found: a 
smaller size of the radius of the rare earth is such to 
decrease the Ni-O-Mn angles increasing the octahedral 
distortions, which in turn reduces the effective interac- 
tions between Ni and Mn ions. Following this procedure, 
the total energies of different magnetic structures as con- 
sidered above for R=La reveal a transition from a FM 
(for R=La,Sm) to the E*-type magnetic structure (for 
smaller rare-earth ions R=Y) [3]. RNMO with R=Sm, 
Y are found to be insulators with electronic structures 
very close to that of LNMO. We show in Fig. Q]b) the 
spin resolved DOS of YNMO in the E*-type magnetic 
structure. Small differences are found in the values of the 
band gaps and magnetic moments at Ni and Mn sites be- 
tween the considered magnetic structures in all RNMO 
studied in this work. The magnetism is governed by the 
superexchange interactions due to the Hund's rule and 



energy gain allows virtual hopping of parallel spins and 
forbids anti-aligned spins between half-filled Ni-e g and 
empty Mn-e 9 orbitals as expected in an extended Kugel- 
Khomskii model 0, E3] • The trend of the exchange pa- 
rameters with changing R is reported in Table HI [18J In- 
plane exchange interactions between Mn sites are negli- 
gible compared to Ni sites. The ratio between the anti- 
ferromagnetic next-nearest-neighbour (NNN), Jj^jvfLjvi) 
and FM nearest-neighbour (NN) Jj^jvi-Mn} interactions 
increases and explains the stabilization of E*-type mag- 
netic structure: a frustrated T~T~4~4 Ising-like spin chain 
in ab planes [l|, [l^l ■ Magnetic couplings changed by in- 
tercalation of rare earth with smaller ionic radius gives 
the possiblity to achieve a polar state in RNMO. We 
recall that the E-type magnetic structure is known to 
produce large polarization along the short b axis due to 
the noncentrosymmetric collinear spin arrangment which 
breaks the inversion symmetry in RMn03 and RNi03 
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In order to understand the transition from the FM to 
the E*-type magnetic state at a more microscopic level, 
we analyze the following two-band model Hamiltonian 
for the double perovskites: 

a/3 
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Here, 



and 



are annihilation and creation 



operators for electrons with spin a in the e g orbital 
a e {x 2 - y 2 (a),3z 2 - r 2 (b)}. denote the hop- 

ping amplitudes between the two e g orbitals on NN and 
NNN sites. The NN Ni-Mn hopping is given by [13|: 
tf 1 = t\ x = t, t%2 = t'22 — tf 2 = ffi = —t/i/3, t\ 2 = 
t\i = t/y/3, where x and y mark the spatial directions. 
The NNN hoppings are parameterized by t' such that 
^a^ N = t' tad 6i denotes the on-site energy, and, 

guided by the electronic structure shown in Fig[TJ we set 
€i(Mn/Ni) = ±A. Jh denotes the Hund's rule coupling 
between the e g electrons and U' (U) is inter (intra)-band 
Hubbard repulsion with U = U' + 2Jh- Since the Mn 
sites have a half-filled t2 S level, an additional Hund's rule 
coupling between the Mn-t2 9 and the Mn-e 9 is included. 
The <T; denotes the electronic spin operator defined as 
°i = c L<x r L' c «*<7', where 1^ are the Pauli ma- 

trices. A denotes the strength of the electron-lattice cou- 
pling. Qxi, Qzi are the Jahn- Teller (JT) distortions which 
are very weak compared to the breathing mode (BM) 
distortions Qu- The r xi = 1 "' 
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Aba c iba) are the orrj ital pseudospin 
operators for the e g orbitals, and the elastic energy for the 
distortions is included as H e i = (1/2) (Q^ + Q 2 zl + Q^). 
Any deviation from collincar spin arrangements is ne- 
glected. 

We identify the key parameters of the model using in- 
puts from DFT results. An important feature from DFT 
calculations is the increase in the relative strength of the 
NNN respect to the NN interactions as R changes from 
La to Y (cfr Table I). Therefore we focus on this varia- 
tion, keeping the other model parameters fixed to reason- 
able values guided by the DFT results 0. A staggered 
pattern of the BM distortions with Qu = Qb e 1 ( 7r,7r '' r * 
is assumed following the DFT results discussed ear- 
lier. A weak JT distortion is also considered in or- 
der to test the robustness of the results. We analyze 
the model Hamiltonian by employing a Hartree-Fock de- 
coupling of the interaction terms; replacing ni a ^rii a i by 
(n ia ^)n iai + n ia ^(n iai ) - {n ia ^) {m a i) , etc. The result- 
ing model is solved self-consistently by iterative numer- 
ical diagonalization on a two-dimensional Ni-Mn square 
lattice starting from different possible magnetic states, 
and the energies of the converged states are compared to 
find the groundstate. In Fig. [2] we show the total energy 
per lattice site as a function of the ratio t' between the 
NNN and the NN hopping strengths. Panel (a) corre- 
sponds to the absence of JT distortions while weak JT 
distortions are present for the results in panel (b). The 
figure demonstrates that it is indeed possible to drive a 
transition from a ferromagnetic to an E*-type state by 
increasing the NNN hopping parameter. JT distortions 
work in favor of the E*-type state as the t' required for 
the FM to E*-type transition reduces upon increasing Q x 
(See Fig. [2b). 

Using the Berry-phase(BP) approach we evaluate the 
first-principles electric polarization P 25[ for YNMO. 



The purely electronic contribution to polarization, P e ; e , 
can be calculated by imposing the E*-type magnetic 
structure in the centrosymmetric crystal structure with 
P2i/n symmetry, leading to P e ; e =0.75 fiC/cm 2 along 
the b axis. The checkerboard pattern due to Ni 2+ and 
Mn 4+ valence in the monoclinic P2i/n crystal structure 
leads to the BM distortions of the NiC>3 and MJ1O3 oc- 
tahedra, which contain three inequivalent oxygen atoms 
(Oi out of plane and O2, O3 in plane) with correspond- 
ing Ni-O-Mn (Ni-0 2 -Mn and Ni-0 3 -Mn) angles differing 
slightly due to the monoclinic distortion of the unit cell 
(see Fig. [3] a)). The E*-type magnetic order further 
lowers the crystal symmetry to P2i. The magnetic or- 
der couples to the lattice and influences the ionic posi- 
tions and the electronic charge distribution, leading to 
a decrease in the electronic contribution with respect to 
the purely electronic term calculated for the centrosym- 
metric crystal structure. However, the total polariza- 
tion increases with respect to the centrosymmetric crys- 
tal structure due to the ionic contribution: Pynmo= 
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FIG. 2: (Color online) Total energy as a function of the ratio t' 
of the NNN to the NN hopping parameters, (a) in the absence, 
and (b) in the presence, of JT distortions. The E-type state 
(equivalent to E* since the model is 2D) becomes lower in 
energy with increasing t' in both cases. The other model 
parameters are A = 2, U' = 4, Jh = 0.5, A = 1, Q z = 0, 
Qb = 1. (c) Total energy as a function of St for A = 2.0, 
J H = 0.5, Q x = and U' = 4. (d) Energy gap in the DOS 
and the net magnetization as a function of St. 
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FIG. 3: (Color online) In-plane arrangement of the ions and 
their relative distances in the unit cell in P2i/c (a) and P2i 
(b) symmetry. Arrows indicate the magnetic spin ordering in 
the E*-type magnetic structure 



2.50 iiC/cm 2 (P e;e =0.51 fiC/cm 2 , P«,„ ic =1.99 fiC/cm 2 ). 
The exchange energy in the P2i symmetry is minimized: 
(i) by increasing (decreasing) the distance between f-f 
(t-|) along the t-T-|-| s P m chains (see d , d in Fig. 
[3]b)), which are identical in the original P2i/n symme- 
try (see d in Fig. [3]); (ii) by further splitting from a 
magnetic point of view of the oxygen ions into 4 inequiv- 
alent sites (see Ni-0 2 --Mn, Ni-0 2 »-Mn and Ni-0 3 <-Mn, 
Ni-0 3 "-Mn in Fig. Ob)). The inequivalent oxygen atoms 
2 ' , 2 " and 3 ' , 3 " are charge-polarized with local 
dipole moments leading to a net polarization developing 
along the b axis. The (2,2) components of Born effec- 
tive charges (BEC) are consistent with the inequivalency 
of oxygen sites: Z* Q , =-2.87 e, Z* Q , =-2.0 e, Z* Q n =-3.65 

e, Zq n =-2.76 e (for the labelling of the oxygen sites 
see Fig. Within the so called point charge model 
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(PCM) Prow ==1.64 nC/cm 2 , i.e. very close to the one 
calculated by DFT, confirming that the electronic rear- 
rangement of the charge due to the symmetry breaking is 
small. Finally we demonstrate using the model Hamilto- 
nian that an external electric field can induce a magnetic 
transition in this system. The main effect of an external 
electric field is to change the on-site energies of the oxy- 
gen, which enter in deriving the hopping parameters used 
in the model. To leading order the hoppings between Mn- 
d and Ni-d via the O-p orbitals are cx tp d /A p d- Due to 
the inequivalence of bridging oxygens( 2 ' an d 2 " in 
Fig. [3]), an external electric field makes A p d and hence 
the d-d hopping parameters staggered. We use a param- 
eter St as a measure of this modulation, such that the 
hopping in both x and y direction are spatially staggered 
(1 ± St) t^g . Presence of an electric field (modelled as 
parameter St) leads to a change in the magnetic ground- 
state of the system from a FM to the E*-type, as seen 
in Fig. |2~fc. Fig. [2ji shows the discontinuous change in 
the total magnetization and the gap in the DOS across 
this transition. The jump observed in the magnitude of 
the gap is consistent with the DFT values for the two 
magnetic phases (see fig. [I]). 

In conclusion, we have performed density functional 
theory and model calculations to elucidate the possibil- 
ity to find novel multiferroics in rare-earths R^NiMnOe 
double-perovskites. The magnetic state changes from a 
FM (R=La) to an E*-type (R=Y): the later being well 
known to assist in ferroelectricity. Calculations on a two- 
band model combined with inputs from DFT calculations 
describe this magnetic transition. Y2NiMn06 is shown 
to be polar in its magnetic ground state with an intrisic 
polarization c omp arable with other magnetically-driven 
ferroelectrics 2(| 27, 2|| 2{|. We suggest experimental 
investigations to test our theorethical predictions. We 
further show that an external electric field can change 
the magnetic ordering from FM to E*-type leading to a 
flip from a non-polar state with finite magnetization to 
a polar state with zero net magnetic moment: a major 
effect that is of interest for electronic devices. 
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